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We study the effect of primordial magnetic fields (PMFs) on the anisotropics of the cosmic mi- 
crowave background (CMB). We assume the spectrum of PMFs is described by log-normal distri- 
bution which has a characteristic scale, rather than power-law spectrum. This scale is expected 
to reflect the generation mechanisms and our analysis is complementary to previous studies with 
power-law spectrum. We calculate power spectra of energy density and Lorentz force of the log- 
normal PMFs, and then calculate CMB temperature and polarization angular power spectra from 
scalar, vector, and tensor modes of perturbations generated from such PMFs. By comparing these 
spectra with WMAP7, QUaD, CBI, Boomerang, and ACBAR data sets, we find that the current 
CMB data set places the strongest constraint at fc ~ 10~^'^ Mpc~^ with the upper limit -6^3 nG. 
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INTRODUCTION 

Many researchers have studied magnetic fields for wider ranges in the early Universe. Recently, magnetic fields 
with strength ^ 1/i G have been observed in clusters of galaxies, and primordial magnetic fields (PMFs) have been 
studied by many authors to explain their origins. The strength of PMFs is constrained as B ^ 1 nG from cosmological 
observations, such as temperature and polarization anisotropics of cosmic microwave background (CMB) and matter 
power spectra [ll-E9|. 

Many theoretical models have been proposed for generating PMFs of cosmological scales. Generation models of 
PMFs from inflation can create strong fields, whose amplitude is of the order of 1 nG at redshift z ~ 0, depending on 
assumed hypothetical fields and couplings 33-34|- Furthermore, PMFs can be generated by primordial perturbations 
of density fields (sHI . l36j , the Beiermann battery mechanism in primordial supernovae remnants (STj or the Weibel 
instabilities [sl]. These PMFs can evolve into the magnetic fields observed in galaxies and/or galaxy clusters directly 
or through the dynamo process. 

While magnetogenesis during inflation can produce PMFs beyond the horizon scale, PMFs generated by the other 
mechanism are expected to have the characteristic scale because they are based on causal processes. Here it should 
be noted that the coherence length of PMFs can grow after generation due to inverse cascade 3^ i^] , although the 
efficiency is still under study. 

In the previous studies to put constraints on PMFs from CMB anisotropies, the spectrum of PMFs was assumed 
to be power-law shape. Although power-law spectrum is natural for magnetogenesis during inflation, this is not the 
case for the other mechanisms based on causal processes. Furthermore, in previous studies with power-law PMFs, it 
was not so clear which scale of PMFs mainly contributes to the constraints. Since inflationary mechanisms tend to 
generate power-law magnetic fields, the magnetic fields which have a characteristic scale at cosmological scales are 
difficult to be produced and seem somewhat artificial. However, it is still possible that such fields can be constrained 
by observation, and if observed it would give us useful information about generation mechanism of cosmological 
magnetic fields. Hence, in this article, as a toy example we use a log-normal distribution (LND) for the PMF 
spectrum /lndI^; ^m, cm) which has a characteristic scale expressed by ku, and aim to constrain PMFs scale by 
scale instead of the power-law PMF spectrum. These quantities refiect the generation mechanism but here we regard 
them just as parameters. We study features of angular power spectra of the CMB from the PMFs with log-normal 
distribution. Finally we constrain the strength of the PMFs for fixed sets of the parameters. 



POWER SPECTRUM OF PMF 



The electromagnetic tensor has the usual form 
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where Ei and Bi are the electric and magnetic fields. Here we use natural units c 
tensor for electromagnetism is 



h ^ I. The energy- momentum 
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The Maxwell stress tensor, cr*'^ is derived from the space-space components of the electromagnetic energy-momentum 
tensor. 
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where a is the cosmic scale factor. 

We consider the case that the spectrum of the PMFs at large scales had been well established before the CMB 
epoch, after their generation and followed by perhaps a complicated time evolution. In this case we do not need to 
consider the time evolution of magnetic fields in our CMB analysis, e.g. nonliner effects (i^l, of the PMF and we 



consider only the dilution due to the expansion of the Universe. Within the linear approximation [41| we can discard 
the MHD backreaction onto the field itself. The conductivity of the primordial plasma is very large, and magnetic 



3 



fields are frozen in tlie plasma [3]. This is a very good approximation during the epochs of interest here. Furthermore, 
we can neglect the electric field, i.e. E ^ 0, and can decouple the time evolution of the magnetic field from its spatial 
dependence, i.e. B(r, x) — B(x)/a-^ for very large scales. In this way we obtain the following equations, 

r*°[EM] = 7'°'' [EM] = , (5) 

-T*'=[EM] =^^' = ^(2S'S'-^''i?') . (6) 



power spectra of PMFs from the log-normal distribution 

We defined a two point correlation function of the PMFs as follows, 

(B\k)B'y\U)) = (2^)3p[pMF](A;)P'J"(fc)<5(fc- fe') , (7) 

where 

P^i{k) = 5'= (8) 
We follow the convention for the Fourier transform as 

f{k) ^ J (fxexp{ik- x)F{x). (9) 

If PMFs have been generated from inflation the power-law model would be a good representation of the magnetic 
field power spectrum. On the other hand, if the PMFs were created through mechanisms other than inflation, the 
spectrum would have a characteristic scale and may not be described as a simple power-law. 

Furthermore, we need to understand how such PMFs cascade from smaller to larger scales to study effects of the 
PMFs created in the early Universe on the large-scale structures. Actually, we can learn behaviors of PMF cascading 
or inverse cascading by simulations [42j, whose results, however, generally depend on cosmological parameters and 
models. Since such simulations generally take too much time and have a limited dynamical range, it is not efficient 
to estimate distributions of the PMF by some simulations when we study cosmology and astrophysics with the PMF 
quantitatively. 

To avoid these problems, we use in place of the power law a LND defined as, 

, ^ 1 f [In(fc) -ln(fcM)lM / X 

/LND(fc;fcM,crM) = 7^ ^xp < } , (10) 

fcCTMv27r [ ^(Jm J 

where /cm is the characteristic scale depending on the PMF generation model, and ctm is the variance. The variance ctm 
in Ea. (|10p expresses how this distribution is expanded (or concentrated) around the characteristic scale /cm. Therefore, 
the variance ctm niay be related to the cascading of the PMFs. Using Eq. pH]) . the power spectrum of the PMF is 
given by 

-P[PMF](fc) = ^/[LND](fc;fcM,crM)- (11) 

We shall determine the coefficient from the variance of magnetic fields in real space, 

{B\x)BHx))l = Bj (12) 

where d — 27r//cM and 

{B\x)B^*{x')) = J d^k J d^k'exp{-ix- k+ ix' ■ k') {B\k)B^*{U) x w{k)) . (13) 

The window function w{k) is given by the rectangular function as w{k) = 1 for k^ < k < k-^- and w{k) ~ for 
otherwise. Here, 

{ '"fcM+"'JM 
/c+ exp(ln/cM +ncrM) = 10 '"d") , 
/e_ = exp (In/cM - ncTM) = 10 '-(i"' , 
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where n gives the range of integration. From Eq.([T3]), we obtain the amphtude of the log-normal distribution as 

(15) 
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As mentioned above, we can neglect the diffusion of the magnetic field at cosmological scales for the age of the 
Universe. Therefore we can assume that the energy of the LND-PMF is not dissipated; in other words, the integral 
of Eq.(ll) over k from k- to fc+ is time- independent. In this case, fc_ and fc+ should be taken to be ± infinites. In 
practice, we take fc_ and fc+ far enough from the peak of the distribution [? ]. 

Using Eq. ([TU]) and the method presented by Ref. [13, , we obtain power spectra from the LND function for each 
perturbation mode as follows. For the scalar mode, the power spectrum of the energy (pressure) of the PMFs is 



8(27r)4a8 



P[EM:SE]ik;ku,cru) ^ g j dk'k'^ j rfC/LND(|fc'|; ^m, o-m)/m(|^ - ^'|; ^m, o-m) 



{(l+C2)fc2-4fcfc'C + 2fc'^} 



life- ife'P 



the power spectrum of the magnetic tension is 
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and the correlation term between the energy and tension is 
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-P[EM:SET](fc; ^M.^m) ^ 4(-27r)4Q8 J ^k' k' J dC/LND(|fc'|; ^M, CrM)/LND(|fc - fc'li ^M, CTm) 



For the vector mode, the power spectrum of the PMFs is; 
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Finally, for the tensor mode the power spectrum of anisotropic stress of the PMFs is found to be 



-P[EM:T](fc; fcM,gM) ^ iQ{^2TTYa^ J dk'k' J dC/LNodfc'l; ^M, 0-M)/LND(|fc - ^'1; fclM, ^m) 
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Since it is difficult to analytically make convolutions of the log-normal distribution functions, we numerically estimate 
them. 



5 



EVOLUTION EQUATIONS 



In this section, we shall summarize the essential evolution equations for each mode. For details, see Refs. 3, 26| . 
In this paper, we use the metric perturbations of the conformal Newtonian gauge defined by Ref. [UliilU We 
also express a scalar potential which is the gravitational potential in the Newtonian limit by ip, and the other scalar 
potential by (j). The vector and tensor potentials in the metric perturbations are expressed by V and H. The subscripts 
"CDM", "b", "7" and 'V" in the equations in this paper indicate cold dark matter, baryon, photon and neutrino, 
respectively 



Boltzmann equations and Hierarchy 



We can understand the time evolution of temperature perturbations and polarization anisotropics of CMB photons 
from the Boltzmann equations. In this subsection, we introduce the Boltzmann equations as follows (in detail see 
Ref. (441), 
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where jig is a number density of free electrons and ut is the Thomson scattering cross section. Here, an index of AI is 
a kind of perturbation mode, and M — 0, ±1, ±, 2 show scalar, vector and tensor modes respectively. The third term 
s\^^^ shows effects of gravitational potentials and scattering by other matters. We can write down nonzero terms of 
iSl^^"* as follows: 
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where tc — arieUT is the differential optical depth of Compton scattering, and V^^ are the anisotropic effects from 
Compton scattering and polarizations given by 



V^^ = — 
10 
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Here, f l^^-* and BY^''' are E and B mode polarizations, respectively, whose evolutions are given by 
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where 
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Scalar Mode 



From Refs. 16, 2^, 43-43, we obtain the same form for the evolutfon equatfons of photons and baryons as in previous 
work 43|-|47|, by considering the Compton interaction between baryons and photons, 



k^^ + 3Hi^ + H4,) = 4^Ga2 {%m:S] (fc, r) - <5ptot } 
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(29) 
(30) 
(31) 
(32) 

(33) 

(34) 

(35) 

(36) 
(37) 

(38) 



where in the second term on the right-hand side of Eq. ([55t is the shear stress of the photons, and L[em:S](^i ''') 
is the Lorentz force. Here 

E[^M:S] (^' = -P[EM:SE] (fc; fcpMF] , cr) (39) 
-^ffiMiS] (^' = -P[EM:SE] (fc; fc[PMF] 70-) + -P[EM:ST] (^i fc[PMF] , ) - ^[EMiSET] (^i fc[PMF] , <^ ) (40) 

2 4 2 

■^[BM:S](^''^) = q-P[EM:SE](fc; fc[PMF],0-) + -P[EM:ST] (^i ^[PMF] , Cr) ^ o -^[EMtSET] (^; ^[PMF] , O") ■ (41) 



Vector Mode 

We can obtain the evolution of the vector potential V{t, k) influenced by a stochastic PMF as [i^, IS] 

k(v + 2-v]= -%TTa?G [2n[EM:V] (k, t) + p^^tt^ + p^Tr^] (42) 



where the dot denotes a confornial time derivative, while pi and tt^ are the pressure and the anisotropic stress of the 
photons {i = 7) and neutrinos (i = v)- Since the vector anisotropic stress of the baryonic plasma fluid is negligible 
generally, it is omitted. Following [ij, |4^, the Euler equations affected by the PMF for the neutrino, photon and 
baryon velocities, v'l^\ and v^'^ are written as 

i'^'-V^-k[^Q^^i^ , (43) 
^(^)-T> + r.(.(v)-«r)) = -fc(^^Qg^ , (44) 

a it 

^^ntEM:V](k,r) ^^^^ 
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where R = (3/4)(pb/p^). Here, the vector Lorentz force is given by 

n^EMiV] (k, t) = P[EM:V] {k; fc[PMF] , ct) . (46) 
For the photons v'^^ — Q^i \ while q[^] and q''^2 quadrupole moments of the neutrino and photon angular 



distributions, respectively. These quantities are proportional to the anisotropic stress tensors. Equations (|43p -(j45 
denote the vector equations of motion for the cosmic fluid, which arise from the conservation of energy-momentum. 

Tensor mode 



The tensor mode H is given by [4J-|47| 



n + 2^H + k^H = 8TTGa^ (n[EM:T]+^''P) , (47) 
where Hi ' is the anisotropic stress for neutrinos. Here, 

n^EMiT] O^^ = ^[EM:T] {k; fc[PMF] , Cr) ■ (48) 
Initial conditions 

We adopt the " compensated Magnetic Modes" as our initial conditions [2^ |4§| . Since these initial conditions have 
been explained in detail in Appendix B of Ref.[2^ and Sec. 4 Ref.[4^, and the list of solutions is too long, we give a 
brief introduction to these in this article as follows. At first a photon fluid is tightly coupled with a baryon fluid in 
the very early Universe. Therefore we assume the Thomson scattering terms are negligible. In this case we can set 
krc <^ 1 and tc/t <^ 1. We assume that the baryon has no pressure and the baryon-photon fluid is representable as 
an ideal fluid and neglect the anisotropic pressure perturbations of the fluid. 

RESULTS AND DISCUSSIONS 

In this section we will discuss dependencies of power spectra from the LND-PMF[Eas. p^ - (pil)) ] on the characteristic 
scale fcjvi and variance ctm, and show results of t emp erature and polarization anisotropics of the CMB generated from 
the LND-PMF with a modified CAMB code 0. We wiU, also, discuss how the strength of the LND-PMF is 
constrained by the CMB observations. 

Effects of the LND-PMF on the CMB 

Figures [D and [2] show the four modes of the CMB (TT, EE, BB, and TE modes) from the LND-PMF. Since 
qualitative features of the scalar, vector, and tensor modes on the CMB are almost the same as previous studies 
22 . 23I 25I 28 1, in this article, we will briefly explain these features of the CMB with the PMF as follows. 



In general, the vector mode dominates all the temperature and polarization auto- and cross-correlation anisotropics 
at £ > 500, while the scalar mode dominates at smaller £, except for the BB mode. For BB mode of the CMB, the 
vector mode dominates the power spectrum for all i. 

To understand the feature of temperature fluctuations and polarization anisotropics of the CMB from the LND-PMF 
in FiglU we shafl go back to the definitions Eqs. lfTOl) . (fTB]) — ([20|) and show Figl3]which depicts P{k) and P{k) x k^. 
From the definition of Eq. (fTO|l . we can understand that the peak amplitudes of the LND-PMF power spectra of each 
ku for the CMB spectra are along the scale invariant amplitude [see panel (c) of Fig[3]. Therefore, the superposition 
of the spectra of CMB anisotropy from different corresponds to the spectra from the scale invariant magnetic 
fields. Indeed in the top-center panel of Fig[Tl we find that the points of the CMB power spectra on the same ctm, 
which are corresponding to the peak of the P{k) x k^, are distributed along the scale invariant curve. We should note 
that the peak position of the LND-PMF source P{k) x k^ is different from one of the eventual power spectra of the 
CMB, because the final spectrum should be obtained by multiplying the LND-PMF power spectrum by the linear 
transfer functions which include various physical effects including the PMF. 
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Considering the mathematical definition of Ea. (|TO)) as mentioned above, the peak scale of the LND-PMF power 
spectrum is determined by ku- We find, however, that the peak positions of the CMB spectra from the PMF power 
spectrum are characterized not only by /cm but also by ctm- The reason is as follows. The LND-PMF power spectrum 
P{k) becomes more broad with increasing ctm as shown in panel (b) of Fig. 3. The amplitudes of the CMB power 
spectra are determined by k^P{k) which is plotted in the panel (d) of Fig. 3. The peak position in terms of the 
combination k^P{k) is affected also by ctm, and hence the positions of the peaks in CMB power spectra are also 
affected by au- 



Considering the scale dependencies of the scalar and vector modes [17|, [18|, |26|, |27[, in the case of the smaller 
/cm and/or smaller cm the scalar mode determines the maximum amplitude of the power spectra of the CMB and 
dominates the power spectra at smaller angular moments £ < 1000. Meanwhile, in the case of the larger ku and/or 
larger ctm, the vector mode determines the maximum amplitude and dominates at larger £. In this case, since the 
peak of the total power spectrum from the LND-PMF lies among 1000 < £ < 2000 (e.g., for ctm = 1-0, fc[LND] ^ 1.0^^, 
see Fig. 4), we can expect that the observation data of the CMB among 1000 < £ < 2000 better constrain LND-PMF 
parameters. 



How are LND-PMF parameters constrained by the CMB? 

In the above sections, we have explained the behaviors of the power spectrum of the PMF from LND distributions 
and the effects on the CMB. In this subsection, considering these features, we discuss how LND-PMF parameters, 
km, Cm and Slnd are constrained by the CMB observations. As mentioned in the above sections, the peak positions 
of the spectra of PMF electromagnetic tensors mainly depend on the characteristic wavenumber kM and um, the 
amplitudes mainly depend on ku and Slnd, and the features the peak width and the damping scale depend on the 
CTM. Figures m and [5] show the CMB (TT, EE, BB and TE) angular power spectra from primary (standard adiabatic 
mode without the PMF) and from the LND-PMF along with observations (WMAP ACBAR [ill, CBI |;53] and 



QUaD [5J|). Considering the above mentioned facts, if ku or ctm is larger, the PMF amplitude is constrained by the 
CMB observation data at smaller scales. On the other hand, if or <tm is smaller, the PMF amplitude is constrained 
by the CMB at larger scales. The strength of the PMF Blnd monotonically increases or decreases the amplitude 
of CMB power spectra from the LND-PMF. These features may suggest a significant anticorrelation between the 
characteristic wavenumber /cm and the variance ctm- Meanwhile, the slopes of the spectrum at small scale are mainly 
characterized by the variance ctm- 

To research how the strength of the LND-PMF is constrained by the CMB observations on each characteristic 
scale /cm, we perform a Markov chain Monte Carlo analysis with the CMB observations and obtain the constraint 
on the -Blnd-^m plane. Since our purpose is qualitatively understanding how the primordial magnetic fields are 
constrained scale by scale the variance parameter of the LND-PMF is fixed at au = 1 for simplicity and we also fix 



the cosmological parameters to those of the best-fit flat ACDM model as given in the WMAP 7yr analysis [51 1. 

Figure H shows the constrained strengths of the LND-PMF from the CMB observations (WMAP 7th, ACBAR, 
CBI, Boomerang, and QuaD) for each scale. From this figure, the LND-PMFs at 10"^ Mpc < k < 10~^ Mpc are 
most tightly constrained. 



SUMMARY 



It is natural to assume that the power spectrum of the PMF from inflation is the power law. However, when we 
consider PMF generated by causal mechanisms, it would be appropriate to assume a spectrum with a characteristic 
scale. In this paper, we assumed the log-normal distribution as the PMF spectrum and studied the effects of the 
PMFs on the temperature and polarization anisotropics of the CMB. We have revealed the features of the LND-PMF 
effects on the CMB as follows: (1) For larger and/or ctm, the CMB spectra of TT, TE and EE from the LND-PMF 
are dominated by the vector mode. Meanwhile, in the opposite case, these spectra are dominated by the scalar mode. 
(2) The CMB spectrum of the BB mode for all scales and other modes for smaller scales is dominated by the vector 
modes. Because three parameters which characterize the LND-PMF affect the CMB power spectra differently at 
small and large scales we expect that tight constraints can be placed on these parameters without degeneracy. We 
found that the LND-PMF which generates CMB anisotropics among 1000 < £ < 2000 is most effectively constrained 
by the current CMB data sets. For example, Blnd for 10~^Mpc~^ < k/h < lO^^Mpc^^ and ctm = 1-0 is limited 
most strongly as shown in FigjSl In the near future, the tighter constraints on -Blnd at ku/h > lO^^Mpc"^ will be 
expected from the observations, such as the Planck, QUIET, and PolarBear missions. 
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FIG. 1. CMB spectra from the LND-PMF at ctm = 1.0 and Blnd =10nG. TT, EE, BB and TE in each panel indicate the 
temperature auto-correlation, the E-mode auto-correlation, the B-mode auto-correlation and the temperature E-mode cross 
correlation, respectively Curves in all panels are theoretical lines as indicated in the legend box at the left bottom. Curves of 
all TE modes are plotted in the absolute value. 
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FIG. 6. Constraint on the strengths of LND-PMF for 10"^ < km < 10"\ The bold curve is the 2 a upper limits of Blnd [nG]. 
We fix the standard cosmological parameters and use the best-fitted value from WMAP 7th + tensor mode[Fu]. 



